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We consider a gravitating spherically symmetric configuration consisting of a scalar field non- 
minimally coupled to ordinary matter in the form of a perfect fluid. For this system we find static, 
f"^ , regular, asymptotically flat solutions for both relativistic and non-relativistic cases, ft is shown that 

£S) ' the presence of the non-minimal interaction leads to substantial changes both in the radial matter 

distribution of the star and in the star's total mass. A simple stability test indicates that, for the 
choice of parameters used in the paper, the solutions are unstable. 
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■ I. INTRODUCTION 

<N . 

In cosmology scalar fields play a key role in models of both the inflation era of the early Universe [l| and the current 
accelerated expansion 0. Scalar fields are also widely used to study smaller scale objects such as the effect of a dark 
jy^ ■ matter, scalar field on the structure of galaxies Q • Moving to still smaller scales there have been studies of boson 
stars [4j - compact (usually spherically symmetric) configurations formed from gravitating scalar field(s). The sizes 
of such configurations varies from the microscopic ( "gravitational atom" ) up to sizes corresponding to massive black 
holes such as those presumed to exist at the center of many galaxies including our own. In this connection, it is 
quite natural to investigate the role that such scalar fields may have in the processes of formation and evolution of 
£3 ' galaxies and more compact astrophysical structure - stars and their clusters. As applied to the stars, it is possible 
to imagine a situation where scalar fields can exist inside usual stars consisting of ordinary matter. The existence of 
such scalar fields would undoubtedly have an impact on the star's inner structure. One such example is the model 
of a star consisting of both ordinary polytropic matter and a ghost scalar field Q . The presence of the field leads 
to the appearance of a tunnel at the center of such stellar configurations which provides new geometric and physical 
properties to such objects. Another example of such stellar objects involves fermion fields interacting with a real 
scalar field $ not only gravitationally but also through the Yukawa coupling of the type — fipipQ, where / is a Yukawa 
' type coupling constant Jh- Further examples of mixed configurations consisting of both boson and fermion fields can 
be found in references 

The investigation in the present work is along the direction of the above works - we consider a scalar field interacting 
\ non-minimally with ordinary matter in the form of a perfect polytropic fluid. An interaction similar to this was used 
y—i • to describe the evolution of dark energy within the framework of chameleon cosmologies [To|, fill . In chameleon 
t-H ' cosmologies the properties of the scalar (chameleon) field depend strongly on the environment [12[ into which it is 
embedded. For example, the mass of the scalar field can change depending on the background environment [l2j . As 
• i— i ' will be shown below, the presence of such an interaction between the polytropic fluid and the scalar field results in a 
, substantial change of the inner structure of the polytropic star. This change is related to the fact that the behavior 
of the scalar field depends the "environment" (i.e. the fluid) into which it is embedded as is the case with chameleon 
cosmology models. Thus we call these configurations chameleon stars. 

In building our model of a chameleon star we will use a real scalar field. In the bulk of the literature complex 
(charged) scalar fields are more often used in constructing models of boson stars. The charge of the complex scalar 
field gencrically leads to repulsion which counteracts the attraction of gravity and thus gives a physical basis for the 
existence of these boson star solutions. Real scalar fields have received less attention because they do not carry a 
charge, and constructing regular stable solutions is more difficult task in this case. The known static solutions for 
real scalar fields are either singular [l]| [HI, or use phantom scalar fields 1^17|. Singular solutions with trivial 



topology and non-phantom fields have been found for both massless fields [13( and for self-interacting fields [14 
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The latter ordinary scalar field solutions are asymptotically flat but the issue of their stability is not quite clear (for 
further discussion see reference [H|). In [l9| regular, static solutions in the presence of the cosmological constant 
were found. In this case the space-time asymptotically approaches de Sitte r sp ace-time. The linear stability analysis 
of these solutions indicate that they are unstable. In contrast the paper [20| found regular, static solutions which 
asymptotically approached anti-de-Sitter space-time, which were stable under a linear stability analysis. In the present 
paper the real scalar field non-minimally interacts with the polytropic fluid. As will be shown below, this leads to the 
existence of regular, asymptotically flat solutions. 

In the absence of gravity, Derrick's theorem [2l|, [22| forbids the existence of non-trivial, stable, regular static D > 3 
dimensional solutions from only scalar fields. This theorem assumes that the potential energy of the system is some 
non-negative function, V(ipi) > 0, vanishing only at its absolute minima. A partial extension of Derrick's theorem to 
the case of general relativity has been studied in [2^, [24[ where it was shown that the existence of asymptotically flat 
particle-like solutions with a regular center and normal (non-phantom) scalar fields is also impossible if V((p) > 0. 
The distinction between the non-gravitational form of Derrick's theorem [2l|, H3| and the gravitational variant 0, [24| 
is that the latter does not touch on the stability of the solutions. 

If one allows other fields besides just scalar fields (e.g. vector or fermion fields) then there are static solutions when 
the spatial dimensionality is three or more. In the absence of gravity, one has the 't Hooft-Polyakov monopole which 
is a finite energy, stable solution for Yang-Mills vector fields plus scalar fields [ID, 26| . In the presence of gravity plus 
a real scalar field there are boson star solutions (for a review see which however have one or more of the bad 
features listed above (e.g. non-trivial topology, not asymptotically flat, not stable) 

In this paper we study the gravitating system of a real scalar field plus a perfect fluid and show that for this system 
there are regular, three-dimensional, static solutions. The most important feature of this system, aside from the 
gravitational interaction, is the presence of the non-minimal coupling between the scalar field and the perfect fluid. It 
is this interaction which allows us to avoid the restrictions of the gravitational version of Derrick's theorem [23|, [24j , 
and find regular static solutions even if the potential energy of the scalar field satisfies V(ip) > 0. We take a polytropic 
equation of state for the perfect fluid. Stellar models using polytropic fluids were investi gate d in detail both within 
the framework of Newtonian gravitational theory [27| and for strong gravitational fields |28j ] . Both of these studies 
hinted at the existence of finite size, regular solutions. These solutions were successfully used to give a description 
of both non-relativistic stars and stars where relativistic effects were important e.g. neutron stars. As will be shown 
below, the inclusion of the non-minimal interaction between the polytropic fluid and the scalar field leads to regular 
solutions which, however, are generally different from usual polytropic stars. 

The paper is organized as follows: In section [TT] the general equations describing a static configuration consisting 
of a real scalar field coupled to a perfect polytropic fluid are derived. In section HTT1 these equations are written for a 
particular case when the potential energy of the scalar field is chosen to have a quadratic mass term and a quartic 
self-interaction term. For this potential energy, in section HVl we give the results of the numerical calculations for this 
potential for both relativistic and non-relativistic cases, and discuss the issue of the stability of the solutions obtained. 
Next, in section [V] we present a simple analytical solution for the non-relativistic case with the scalar field taken to 
be massless and for a special choice of the coupling function. Finally, in section IVII we summarize the main results 
and give some speculations about the physical applications of these chameleon star configurations to neutron stars 
and living stars (i.e. stars still on the main sequence). 



II. DERIVATION OF EQUATIONS FOR A STATIC CONFIGURATION 



As discussed in the introduction we consider a gravitating system of a real scalar field coupled to a perfect fluid. 
The Lagrangian for this system is 



L — —- 



R 



16ttG 



V(ip) + f(<p)L r 



(1) 



Here <p is the real scalar field with the potential V(ip); L m is the Lagrangian of the perfect isotropic fluid i.e. a fluid 
with only one radial pressure; /((/?) is some function describing the non-minimal interaction between the fluid and the 
scalar field. The case / = 1 corresponds to the absence of the non-minimal coupling, but even in this case the two 
sources are still coupled via gravity. 

We choose the Lagrangian for the isentropic perfect fluid to have the form L m = p 0, [30[ . Using this Lagrangian, 
the corresponding energy- momentum tensor is (details are given in Appendix [XJ 



T? = f [(p + P )u t u k - 5fp] + dup&tp - Si 



(2) 
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where p and p are the density and the pressure of the fluid, u l is the four-velocity (here and throughout the paper we 
set c = 1). We take the static metric of the form 



ds 2 = e^dt 2 - e x ^dr 2 - rW 



(3) 



where dil 2 is the metric on the unit 2-sphere. The (§) and (\) components of the Einstein equations for the metric 
([3]) and the energy-momentum tensor ^ are 



G°o = -e- 
G\ = -e 



r 2 r 
1_ z/ 

r 2 r 



8ttGT ° = 8ttG 
SnGTl = 8vrG 



fp+-e- x <p' 2 + V(<p) 

-/ p -i e -v 2 +n^ 



(4) 
(5) 



The equation for the scalar field (p coming from the Lagrangian (TTJ) is 



1 9_ 

i—g dx l 



-99 



ik 



dip 
dx k 



dV df_ 

dip m dp 



Using this field equation above with the perfect fluid L m = p, and the metric ([3]), gives the following scalar field 
equation 



2 - + \{ v '-X') 
r 2 



p' = e> 



dV_ dj_ 

dip dip J 



(6) 



where the prime denotes differentiation with respect to r. The Einstein field equations are not all independent because 
of the relation T* k = 0. The i = 1 component of this equation has the form 



Or + 2 [l1 lo) 



V + - 



Tt--(Ti + T$) 



(7) 



Taking into account the expressions 



T 2 =Ti = -fp+-e- x ip' 2 + V(ip) 



and Tq, Tl from (j4|), ([5]) and using ([6]), allows us to write (J7|) as 



dp 
dr 



1 , , dv 
2 ip + P) d-r- 



(8) 



The matter Lagrangian used here, L m = p, is not the only possibility. Other variants can be found in [31| . However, 
the choice L m = p has the simplifying feature © does not contain additional terms involving the coupling function 

/(¥>)■ 

For a polytropic equation of state one has 



p = kp 1 ', 

where k, 7 are constants. Now one can introduce the new variable 9 defined as [3j 

P = Pc6 n . 



(9) 



(10) 



Here p c is the central density, and the constant n, the polytropic index, is related to 7 via n = 1/(7—1). Putting 
these definitions together gives equation © in the form 



p = k p 1 = kp 



1+1/7 



kp 



1+1/n/m+i 



Using (fTTj) in equation ((5) leads to 



2a (n + 1 — = - {1 + aO — , 

ctr ar 



(11) 



(12) 
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with <7 = kp~J n = p c / p c and p c is the pressure of the fluid at the center of the configuration. This equation may be 
integrated to give e v in terms of 0: 



l + o- 

1 + C70 



2(n+l) 



(13) 



where e v ° is the value of e v at the center of the configuration where 0=1. The integration constant v c , corresponds 
to the value of v at the center of the configuration. It is determined by requiring e v = 1 at infinity that i.e. that the 
space-time is asymptotically flat. 

The gravitating system of a real scalar field interacting with a perfect fluid is characterized by three unknown 
functions - A, and ip. These three functions are determined by the three equations (|4|). © and ©, and also by the 
relation (|13[) . We now rewrite these equations by the introduction of a new function u(r) [28j 



u(r) 



2GM 



(1 



1 



2GMu 



(14) 



Here M is the mass of the configuration within the range < r < r^, where r& is the boundary of the fluid where 
= 0. Using this function, equation (|4]) becomes 



M— = 4irr 2 
dr 



2GMu\ 



tp' 2 + V 



(15) 



From (|14j) one can define M(r) = Mu(r) which can be interpreted as the total mass of the configuration in the range 
[0, r\. This mass has contributions from both the fluid and the scalar field, within a sphere of coordinate radius r. To 
avoid a singularity in M(r) at the origin, one has to put u(0) = [28[. This corresponds to the fact that the mass at 
the origin is equal to zero i.e. M(0) = 0. 

In anticipation of analyzing the system of equations — (I3J) i (EJ) and ©) - numerically, we introduce the following 
dimcnsionlcss variables 



£ = Ar, v(0 



A 3 M 



u(r), m 



4ttG 



1/2 



<p(r), where A 



4irGp c 



1/2 



.*(n+l)J rv " " " [(n + ^kpV" 

A has the dimensions of an inverse length. With this one can rewrite equations Q and ([5]) in the form 

2 ~1 



dr 



di=Z \ fen+ 2 



l-2<7(n + l)- 



2a(n+l)v 



dO 



l + a6l d£ 



f9 n (1 - 0-0) + 2V 



1 dv 

¥di 



(16) 

(17) 
(18) 



where V = V/p c is the dimensionless potential energy of the field. 
Next, using (|12|). one can rewrite equation §Q as follows: 



de 



a(n + 1) 

1 + CT0 



1+0-0 1 f dv 



1 - 2a(n + l) 



de 



7^ af) n+ldf\ 



di 



(19) 



Thus the static configuration under consideration is described by the three equations (fTT)) - (fll?)) . 



III. CONFIGURATION WITH A MASS AND A QUARTIC SELF-INTERACTION TERM 

In this section we will show that there are non-singular, finite- mass solutions of equations ()17[) - (|19p . First we 
specify the boundary conditions. Using the above dimensionless variables we "normalize" to unity at the center of 
the configuration £ = 



0o = 0(0) = 1. 



(20) 
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From equation (fl7|) one can show that v —> like £ 3 as £ — > 0. Combining this with equation (|18|) one in turn finds 
d9/d£ — > as £ — > 0. Bearing in mind that we are looking for regular solutions, we define the boundary conditions in 
the vicinity of £ « as 



2* ' 



v 3 e 



(21) 



where cf>o corresponds to the initial value of the scalar field (j), the parameters 82 , v 3 are arbitrary, and the value of the 
coefficient 4>2 is defined from equation (|T9)) as 



1 



iii+i 



4f 
#7o 



The index denotes that the values of the functions are taken at £ = 0. 

Using the boundary conditions ([2~Tj) , we proceed to solve the system l|17 p -(|19 p numerically. The behavior of the 
solution will depend both on the parameters of the polytropic fluid n, u and the form of the potential energy of the 
scalar field V, and the coupling function /. One of the simplest and most commonly used choices for the potential 
energy is that of a scalar field with mass m and a quartic self-interaction term (n > 0) 



V 



- ™2 2 

-m ip 



-Kif 



It was pointed out in [14[ that the system with such potential has only singular static solutions. Below we show that 
an inclusion of the non-minimal interaction between the scalar field and the polytropic fluid allows one to obtain 
regular static solutions. As an example, let us choose the coupling function / in dimensionless form as follows 



Then, using the dimensionless variables introduced in the previous section, the potential can be rewritten as 



1 , 

where the new dimensionless constants are given as 

m 2 o~{n + 1) 



A 



K 
Pc 



(22) 



(23) 



AttG 



Using the above potential V and the function /, equations ([T7| - ([T9f take the form 



2 1- 2a(n + l)v/£ d9 



- v 2 9 n (l-ad) +/i 2 2 + -A^ 4 - 



d£2 




1 dv 



(24) 
(25) 



l-2o-(n+l)~ 



(26) 



The parameter /j, can be absorbed by introducing the rescalings x = v = fj,v, f) = /3//i 2 , A = A//i 2 . Bearing 
this in mind, we assume fi = 1 in further calculations. One can see from these equations that the presence of the 
interaction between the scalar field and the fluid is defined by the parameter /3. In the absence of the fluid, this 
system only has singular solutions (l4j . The inclusion of the fluid changes the situation since the presence of the term 
(3a9 n+1 on the right hand side of ([26)) for the scalar field corresponds to an effective mass term [fi 2 — (3a9 n+1 ) whose 
sign depends both on the behavior of the fluid density 9 and the values of the parameters (3 and a. We note that 
the present system of gravitating real scalar field plus perfect fluid turns out to be similar to the gravitating complex 
scalar scalar field system considered in [33[ . The gravitating complex scalar field studied in [33| had a potential of the 
type (|23p and it gave regular, stationary solutions. Thus this is already a hint that we can expect regular, stationary 
solutions for the present system. 
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IV. NUMERICAL RESULTS 
A. Relativistic case 

The results of the numerical calculations for the system (|24[) - ([26|) with the boundary conditions (|21|) are presented 
in tables U and |TTJ The solutions were started near the origin (i.e. near £ w 0) and solved out to a point £ = £1 where 
the function 9 becomes zero. From equations © and (fTU)) this is where the fluid vanishes and it is this point, £ = £i, 
that we define to be the surface of the star. Beyond the point £ = £i we continued the numerical solutions with 
only the gravitational and scalar fields while the fluid was set to zero. The interior and exterior solutions were then 
connected to one another. Details of this are given below. Previous numerical studies with a fluid having a polytropic 
equation of state [28| found regular, relativistic star-like configurations. These star-like solutions of [28| where found 
for the values of the parameters < a < 0.75 and 1 < n < 3. Since our system has two additional parameters 
(A,/3) we restrict ourselves to examine just two sets of the parameters (i.e. a = 0.2, n = 1.0 and a = 0.2, n = 1.5) at 
different values of A and (3 when looking for regular solutions. Below we show that these two sets of solutions differ 
considerably in the behavior of their physical characteristics. 

The parameters of the system under consideration for the chosen values of a and n are given in tables HI and HT1 The 
procedure for finding these parameters is the following: given the values of a and n, we seek eigenvalues (i.e. values 
of the parameters A and /?) for which the function 8 goes to zero at some finite value of £ = £i which as mentioned 
above we take to correspond to the surface of the star. Then since M(r) = Mu(r) (where M is the total mass of the 
star between r = and its surface at r = r^) it is required that u{rjj) = 1. Using this we can evaluate the function 
v(£) from (mi) at £ = £ x 

= (27) 

This quantity defines the total mass M of the configuration through the parameters of the fluid a, n and p c which 
determine the parameter A. 

The value £j of the coordinate £, corresponding to the boundary of the configuration, does not represent the radius 
of the star as measured by a distant observer. To define this radius it is necessary to make a coordinate transformation 
to a new dimensionless variable £ which is defined as follows 

e X/2 dt 



£= / e [l-2a(n+lM£)/£r 1/2 d£. (28) 



or taking into account (114)) and ([16 



o 

Then the observable radius of the configuration R in dimensional variables is defined as R = £i /A in accordance with 
the data presented in tables [T] and [TTJ 

Also the ratio of the central density p c (of the fluid only) to the average density p (of both the fluid and the scalar 
field) is presented in the tables. We define the average density as [28[ 

M 3MA 3 



3 ■ 



r (4/3)7ri? 3 4^ 

Using the expression for i;(£) from ([27]) . one can express the quantity MA 3 /4ir in terms of the central density p c and 
the boundary value u(£i) of the mass function. Then we obtain the relation 

p c _ £? 



p 3w(&) 

Finally, in the last column of both tables the gravitational potential energy expressed in units of M is shown (see 
[28| for more details): 

n 1 i /* r °£ 2 d£ 



M ~ u(£i) J Q [i - 2(7(77. + 1M0/£] 1/2 ' 
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TABLE I: The parameters of the relativistic configurations with the constant central value of the scalar field <f>o = 0.45 



A 


P 


6 


fi 




pc/p 


-Q./M 


a = 0.2, n = 1.5 


1500 


10200 


0.37 


0.4017 


0269 


0.6273 


0.1344 


1250 


8500 


0.41 


0.4450 


0295 


0.7786 


0.1344 


1000 


6900 


0.47 


5092 


0396 


1.0609 


0.1341 


750 


5200 


0.54 


0.5851 


0.0376 


1.3972 


0.1341 


500 


3500 


0.66 


0.7150 


0.0458 


2 0939 


0.1340 


250 


1800 


0.95 


1.0277 


0.0637 


4.4834 


0.1337 


150 


1130 


1.25 


1.3494 


0.0801 


8.1296 


0.1331 


100 


790 


1.58 


1.7012 


0955 


13.7681 


0.1326 


75 


620 


1.83 


1.9678 


0.1075 


19.0029 


0.1322 


60 


525 


2.13 


2.2830 


0.1160 


27.7771 


0.1316 


41 


395 


2.70 


2.8820 


0.1330 


49.3308 


0.1308 


30 


325 


3.35 


3.5573 


0.1449 


86.4977 


0.1299 


20 


260.5 


4.70 


4.9472 


0.1594 


217.0848 


0.1287 


10 


195.24 


9.00 


9.3259 


0.1798 


1351.2011 


0.1269 


5 


162.1683 


21.00 


21.4257 


0.1983 


15568.1073 


0.1229 


3 


148.8345216 


30.00 


30.4723 


0.2102 


42819.2173 


0.1192 


a = 0.2, n = 1.0 


1500.0 


8700 


0.210 


0.2315 


0.0313 


0.0985 


0.1229 


1250.0 


7250 


0.230 


0.2535 


0.0343 


0.1181 


0.1229 


1000.0 


5800 


0.257 


0.2833 


0.0384 


0.1473 


0.1230 


500.0 


2910 


0.362 


0.3992 


0.0543 


0.2909 


0.1231 


350.0 


2050 


0.430 


0.4742 


0.0648 


0.4090 


0.1231 


200.0 


1200 


0.567 


0.6250 


0.0846 


0.7181 


0.1229 


150.0 


920 


0.655 


0.7216 


0.0965 


0.9712 


0.1227 


100.0 


640 


0.800 


0.8805 


0.1152 


1.4810 


0.1222 


50.0 


350 


1.110 


1.2204 


0.1559 


2.9250 


0.1217 


30.0 


238 


1.400 


1.5362 


0.1880 


4.8644 


0.1209 


20.0 


180 


1.650 


1.8084 


0.2160 


6.9318 


0.1203 


10.0 


126 


2.160 


2.3567 


0.2551 


13.1668 


0.1185 


6.8 


109 


2.442 


2.6577 


0.2724 


17.8189 


0.1175 


5.0 


100 


2.670 


2.8992 


0.2825 


22.4626 


0.1167 


3.0 


90 


3.000 


3.2474 


0.2950 


30.5069 


0.1157 


1.0 


80 


3.500 


3.7719 


0.3094 


46.1887 


0.1145 


0.3 


76 


3.670 


3.9511 


0.3169 


52.0002 


0.1141 


0.1 


75 


3.750 


4.0345 


0.3186 


55.1796 


0.1140 



where the expression for the energy density is defined from Using equations (|T^|) . JT 
write the energy density in the following form (in units of p c ) 



and 



-L n 




1 - 2a(n + 1) 



one can 



(29) 



The absolute value of the potential energy H represents the work that would have to be done on the system to diffuse 
its mass to infinity. 

Next, using the data from the tables, in figure [T] we plot f(£i), which is proportional to the total mass of the 
configuration M, on the value of the self-coupling parameter A at some constant central </> = 0.45. In figure [5] wc 
plot v(£i) as a function of the value of central cf>o for different n values and a = 0.2. From figure Q] one can see that, 
while A increases, the masses of the configurations decrease both for n = 1.0 and n = 1.5. Initially, for small A, the 
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TABLE II: The parameters of the relativistic configurations as a function of the central value of the scalar field (j>o for fixed 
observable radius of the star £1 and for given values of a and n. 



00 


A 


P 






Pc/p 


-il/M 






a = 0.2, n = 1.5, £i = 


3.219 










Relativistic star without 


a scalar field 












2.699 


0.9604 


6.8270 


0.2207 


Relativistic star with the scalar field 


0.30 


2101.0 


1642.5 


3.055 


0.1022 


92.9763 


0.1325 


0.35 


617.0 


868.0 


3.040 


0.1186 


78.9467 


0.1323 


0.40 


184.0 


531.0 


3.030 


0.1306 


70.9796 


0.1316 


0.44 


54.5 


385.0 


3.025 


0.1376 


67.0397 


0.1308 


0.45 


34.9 


355.0 


3.025 


0.1396 


66.0795 


0.1304 


0.47 


6.6 


317.0 


3.025 


0.1402 


65.8043 


0.1297 






a = 0.2, n = 1.0, £i = 


2.657 










Relativistic star without 


a scalar field 












2.277 


1.1430 


3.4430 


0.2108 


Relativistic star with the scalar field 


0.27 


510.15 


275.10 


2.424 


0.2986 


15.8983 


0.1262 


0.30 


300.40 


225.90 


2.425 


0.2947 


16.1309 


0.1250 


0.33 


178.90 


190.30 


2.427 


0.2903 


16.4159 


0.1237 


0.35 


124.50 


169.95 


2.430 


0.2889 


16.5532 


0.1228 


0.39 


57.70 


142.00 


2.436 


0.2794 


17.2429 


0.1207 


0.41 


36.00 


130.00 


2.442 


0.2769 


17.5278 


0.1197 


0.43 


19.30 


119.00 


2.442 


0.2739 


17.7255 


0.1185 


0.45 


6.80 


109.00 


2.442 


0.2724 


17.8189 


0.1175 


0.46 


1.90 


105.00 


2.442 


0.2714 


17.8864 


0.1171 



masses are very different for n = 1.0 and n = 1.5; asymptotically, for large A, they approach comparable values. This 
is because for large A the total mass is defined by the scalar field but not by the contribution from the fluid. The 
sizes of the configurations are different for the two cases n = 1.0 (which has a size characterized by £ = 2.657) and 
n = 1.5 (which has a size characterized by £ = 3.219). 

On the other hand if one fixes the size of the star to equal the size of a corresponding relativistic configuration 
without a scalar field, and changes the central value of 0o so as to get regular solutions, it is necessary to choose 
appropriate eigenvalues of A and j3. In this case v (£i) as a function of 4>o is given in figure [2j From this figure one 
sees that as (f>o increases the mass of the configuration can either increase (for the n = 1.5 case) or decreases (for the 
n = 1.0 case). 

It follows from tables U [TT| and figures [TJ [5] that the masses of the configurations, computed for the values of the 
system parameters used here, are considerably smaller than the masses of the stars without a scalar field. The 
maximal dimensionless masses of the configuration with the scalar field for n = 1.0 and at n = 1.5 are 0.32 and 0.21, 
respectively; without a scalar field the masses of the same configurations are 1.14 and 0.96, respectively. Our attempts 
to increase the masses of the system by varying 4>q were not successful. From figure [5] the implication is that the 
masses change only slightly with 4>q. 

Next, in figure [3] the metric functions e A , e v : and the mass distribution M(£/£i) (as a function of the dimensionless 
radius £) are given. In plotting these functions we used the following procedure: 

Inside the Star: The function e v was plotted using (|13[) ; the function e A was plotted from equation (|14p . in terms 
of the dimensionless variables from (fl~6)). Explicitly 

(30) 

From the interior solution parts of figure [3] (i.e. the region < £ < £i) we see the reason for the terminology 
"chameleon star" for the present solutions - the scalar field mimics the behavior of the fluid in the interior region. 
Even in the exterior region where the fluid goes to zero the scalar field is asymptotically going to zero. 



l-2<7(n + l)- 
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0.45 <l>o 



FIG. 1: The dependence of the function v(£i) (which is pro- FIG. 2: The dependence of the function v(£i) on the central 

portional to the total mass M) on the value of the coupling value of 0o at the fixed fi = 3.219 for n = 1.5 and fi = 2.657 

parameter A from (f23|) at the fixed central value <f>o = 0.45 and for n = 1.0, and a = 0.2 for both graphs. The data are taken 

a = 0.2. The data are taken from table [TJ from table ITT1 
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FIG. 3: The density p and the pressure p of the fluid, the scalar field 4>, the metric functions e" , e A and the current mass M 
as functions of the relative invariant radius for n = 1.5 and n = 1.0 respectively. To provide the asymptotic flatness of 
the solutions, i.e. e",e A — > 1 as £ — > oo, the value of the constant v c needed to be chosen as v c ~ —0.958 for n = 1.5 and 
v c w -0.822 for n = 1.0. 



Outside the star. The solution goes to the Schwarzschild solution 



= 1 - 2a(n + 1 



(31) 
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FIG. 4: The distributions of the total energy density in FIG. 5: The distributions of the total energy density in 
units of p c from (|29|l for different the initial values of <po = units of p c from (|29|l for the different initial values of 0o = 
0.30, 0.35, 0.40, 0.45, from top to bottom. 0.45, 0.35, 0.27, from top to bottom. 

The mass function is defined as 

Finally, in figures 2] and [5] we plot the energy density, Tq, from the expression (|2T)|) . The definition of the mass 
function given just above takes into account the energy of the fluid and the scalar field from £ = to £ = £i. Although 
the fluid vanishes at £ = £i the scalar field does not vanish. Thus in principle one should include the contribution to 
the mass function of the scalar field from £ ~ £i to £ —> oo. However from the figures [3][5] one can see that <f> and </>/ 
rapidly go to zero as £ — > oo and thus the mass function is given essentially just by the fluid and scalar field energy 
density between £ = and £ = £i. 

We now turn the the question of connecting the interior region (i.e. < £ < £i where one has both fluid plus scalar 
field as a source) with the exterior region (i.e. £ > £i where one has only a rapidly vanishing scalar field). For this 
purpose, we write down the Einstein equations of © and ([5]), and the scalar field equation of ©, without the fluid 
source i.e. 6 = so that from ^ and (fTT))) we have zero pressure and density. This leads to the following system of 
equations: 

(32) 
(33) 
(34) 



•r 



= 8ttG 

= 8nG 

i > 

(p = e 



-e" V 2 + V(tp) 



V(<p) 



dip ' 



which, by analogy with the transformations made above, can be rewritten in terms of the dimensionlcss variables 
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v (Qi4>(Q an d f(0 as follows 



dv 



(14> 



V 



dv 2a(n+l)e x 




-v 



dv 



dV 



(35) 
(36) 
(37) 



where e A is given by (|30[) . and V is taken from 



Note that this system of equations ([3"2 j) -(|3"I j) (or ([55 |) -(f57 |) ) is 



essentially the same as those obtained in [14| which considered a real scalar field with quartic self-interaction. The 
above system contains the parameter a as a trace of the influence of the fluid on the external solution. The solution 
sought beginning from the surface of the star at £ = £i using, as the boundary conditions, the values of v(£i), 
and obtained from the solution of the equations ([2"4"1) - (|2T)1) for the internal part of the configuration. This allows 

one to determine the value of the integration constant v c from (fT3|) by requiring e u to be equal to unity at infinity, 
providing asymptotical flatness of the space-time. (The values of v c for the examples shown in figures [3] are given 
in the caption.) Thus the complete solution for the configuration under consideration is derived by matching of the 
internal fluid solutions given by equations (f2~4")) - (f2l3)) with the external solutions obtained from the system (|3"5|) - (|3"T1) . 

The system (f3"5")) - (j3"T)) has obvious asymptotically flat solutions in the form: v —> const ps u(£i) with taken 
from (|3"Tj) that corresponds to the fact that, by choosing the eigenvalues of the parameters A and j3 presented in tables 
U and HIl the external scalar field makes negligible contribution to the total mass of the configuration. It allows using 
the solution for the metric functions in the form of (|31[) as a good approximation; v — > const, where the constant is 
made equal to zero by the corresponding choice of v c (see above); <j> — > const e~ M ^/£. The numerical calculations of 
the system ([3"5j) - ([3"T]) confirm this asymptotic behavior. 

From figures [3j [4] and [5] we can draw the following conclusions about the metric functions, e ,e v and the mass 
distribution M(£/£i): 

(1) Asymptotically, as £ — > oo, the space-time becomes flat, i.e. e v , e A — s- 1. To get this asymptotic behavior of the 
metric function e", it is necessary to choose a specific value for the central value, v c (see the caption of figures [3]). 
The asymptotic behavior of e x — > 1 follows from equation ([50]) and the fact that i^)^-^ - > as £ 3 . 

(2) As shown in [28j | . the relativistic configurations without a scalar field are characterized by a greater concentration 
of matter toward the center than in the non-relativistic case. By including a scalar field, we obtain even greater 
concentration of matter toward the center. This can be seen by comparing figures [3] with figures 1 and 2 from the 
paper 28j . The main reason for the increased concentration of mass is the non-minimal coupling between the scalar 
field and the fluid with the choice of the coupling function given by (|22[) . The functional form of this coupling gives 
the required regular solutions, lying in the range < £ < £i, only for the eigenvalues of the parameter ^> 1 (see 
tables HI |TT| . This leads to the fact that near the origin, when £ ss 0, a greater concentration of mass occurs due to 
the presence of the term j3(j) 2 9 n in the expression for the energy density (f2T))) . At the same time both (j> and <f>' give 
relatively small contributions to the energy density compared to the term coming from the non-minimal scalar-fluid 
coupling. This fact will become important, below, when we consider the stability of the solutions. 

(3) The total energy density presented in figures [4] and [5] strongly depends on the polytropic index n. At comparable 
central values of the scalar field 0o the energy density for the case n — 1.5 is several times greater than for the case 
n = 1.0. In turn, both values of n yield a greater energy density at the center of the configurations than in the case 
of the relativistic stars without a scalar field [28j when the energy density was found to be 1. 

(4) Despite the higher concentration of matter at the center of the configurations with the addition of the scalar 
field, their masses are considerably smaller than the masses of the relativistic stars of the same size but without a 
scalar field (see table Hi)) . This occurs as a consequences of the fact that the external region of the stars, with only 
the scalar field as a source, is strongly rarefied due to the rapid vanishing of the scalar field in the region £ > £i. 



B. Non-relativistic case 



In this section we consider the non-relativistic limit of the system (f2"4"l) - ({2"6"]) . The non-relativistic limit corresponds 
to a — > 0, and p -c P- First, we recall some results about this system in the absence of the scalar field [32] ]. With no 
scalar field the system of equations ([24)) - (f2"6"j) reduces to the well known Lane-Emden equation 



1 d 



- 2 d9_- 



(38) 
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This equation has solutions which describe finite size configurations for different values of the parameter n [32| . The 
non-relativistic limit for the case when there a scalar field is obtained by omitting all terms with a in equations 
(|24l) - ([26l) except for the term j3a9 n+1 in the scalar field equation (|26|) . This term must be kept since the product f3a 
is generally a non-zero quantity even for a — > 0. Taking all the above into account we rewrite the system (|24 p - (|26[) in 
the following form: 



,d6_ 

dv 

2 



e 
e 



i 



1 dv 



2 an 



(39) 
(40) 
(41) 



Here, as in the relativistic case the presence of the term f3ad n+1 is important - it leads to a change in sign of the 
effective mass term (/i 2 — f3a8 n+1 ) for certain values of f3. This feature is important for the existence of regular 
solutions. 

Proceeding as in the previous section, we obtain numerical results for the non-relativistic limit which we present in 
table Mil We chose the parameter a to be a = 0.001. Such a small value of a requires that j3 be large enough so that 
the first term on the right hand side of equation (|4"Tj) is negative, which is the necessary condition to have regular 
solutions. In the non-relativistic case, as in the relativistic case, we need to find eigenvalues of two parameters j3 and 
A in order to have regular solutions for a given values of the polytropic index, n, and the central value of the scalar 
field, 4>o- From table Hill we do indeed find that there are values of (3 large enough to make the first term on the right 
hand side of ([4Tj) negative thus yielding regular solutions. 



TABLE III: The parameters of the non-relativistic configurations with a variable central value of the scalar field 4>o- The 
observable radius of a star fi is fixed. The choice of the parameters a and n is indicated in the table. 



4>o 


A 








Po/p 


-il/M 








a = 0.001, n = 1.5, £i = 


= 3.6560 






Non-relativistic star without a scalar field 








3.6538 


2.7141 


5.9907 


0.0000 


Non-relativistic star with the scalar field 


0.039 


1040 




20600 


3.6538 


0.4384 


37.0928 


-0.0009 


0.038 


1850 




21800 


3.6538 


0.4373 


37.1840 


-0.0009 


0.036 


3980 




24500 


3.6538 


0.4354 


37.3411 


-0.0009 


0.034 


7010 




27600 


3.6538 


0.4346 


37.4105 


-0.0009 


0.030 


17575 




35500 


3.6538 


0.4349 


37.3895 


-0.0009 


0.028 


25895 




39500 


3.6538 


0.4462 


36.4388 


-0.0009 


0.027 


31270 




41700 


3.6538 


0.4556 


35.6855 


-0.0009 








a = 0.001, n= 1.0, £i = 


= 3.1437 






Non-relativistic star without a scalar field 








3.1416 


3.1416 


3.2899 


0.0000 


Non-relativistic star with the scalar field 


0.039 


190 




10700 


3.1416 


0.6568 


15.7364 


-0.0008 


0.038 


595 




11000 


3.1416 


0.6687 


15.4556 


-0.0008 


0.036 


1645 




12000 


3.1416 


0.6761 


15.2872 


-0.0008 


0.034 


3070 




13100 


3.1416 


0.6868 


15.0493 


-0.0008 


0.030 


7645 




15700 


3.1416 


0.7166 


14.4227 


-0.0008 


0.025 


20685 




21000 


3.1416 


0.7457 


13.8609 


-0.0008 



Using the data from table [TTTl the dependence of the on t ne central value of the scalar field <f>Q is given in 

figure [HI Figure [7] uses the date from table [TTTl to plot the ratio of the energy density to the central density, T§ / p c as 
a function of the dimensionless, normalized radius. Also from table [TTTl one can see that we have fixed the radius of a 
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FIG. 6: The dependence of the function on the central value 
of </>o at the fixed £1 = 3.6560 for n = 1.5 and £i = 3.1437 for 
n = 1.0, and a = 0.001 for both graphs. The data are taken from 
table ITlTI 



FIG. 7: The plot of the total energy density, T§, in units of 
p c in the non-relativistic case (a = 0.001) for different cen- 
tral values of 4>o from tablc lllll The solid lines correspond to 
the case of n = 1.5, the dashed lines to the case of n = 1.0. 
The top two sets of curves correspond to the configurations 
with the scalar field. For comparison the energy densities 
of the non-relativistic configurations without a scalar field 
Tq = p c 9" are also shown by the two bottom, less peaked 
curves. 



star R = £i to be equal to the size of a non-relativistic star without a scalar field. From the figure we can draw the 
following observations: 

(1) The masses of the stars with a scalar field are considerably smaller than the masses of the stars without a scalar 
field. In the case n = 1.0 the masses are 4-5 times smaller, and in the case n = 1.5 about 6 times smaller. 

(2) The mass is a fairly slowly varying function of </>o, especially when n = 1.5. 

(3) There is a greater concentration of the mass density towards the center of the stellar configurations. This can 
be seen explicitly by comparing with the distribution of the mass density for a star without a scalar field which is 
also presented in figure [7J 

All these results are very similar to those of the relativistic case. 



C. Stability of the solutions 



In this section we discuss the issue of the stability of the regular solutions obtained above. There are two basic 
approaches to studying the stability: (i) The energy approach to the theory of equilibrium for a star (ii) A 
more rigorous dynamical stability approach based on studying the stability of linear and nonlinear time-dependent 
perturbations. In this paper we will use the first approach. Proceeding along the lines of reference (28j . we define the 
total energy E of the system, including the internal and gravitational energies, as 

r R 

E = M = 4tt / T$r 2 dr, (42) 
Jo 

where T§ is defined by (0) and corresponds to the total energy density of the system (recall that c = 1). Next, we 
consider a system consisting of a gas of particles having a rest- mass density p g . Its relativistic energy density is the 
sum of the rest energy p g plus the density of internal energy. For the special case of an adiabatic process which 
assumes the absence of heat flow terms in the energy- momentum tensor ^ , it is possible to obtain a relation between 
the gas density p g and the total mass density p t = Tq or alternatively between p g and 8, tp. 
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To do this we use the first law of thermodynamics which, in our case, takes the form 

dpt + (pt +Pt)^yr = 0, 



(43) 



where pt = ~Tl is the total pressure, and V is the specific volume (this symbol V should not be confused with the 
potential energy used earlier). Since dV/V = —dp g /p g , equation (|4"3")l gives the following relation 



d()g_ 

Pa 



dp t 



Pt+Pt 



(44) 



In general, when p t and pt are functions of 9, ip, this equation cannot be integrated. However, for our configuration, 
the numerical calculations performed in the previous sections show that the main part of the energy is provided by 
the term containing the non-minimal coupling but not by ip' and V(ip). This allows us to neglect the terms containing 
the scalar field kinetic energy and potential as compared with fp and fp in the T§ and T-j 1 components of the energy- 
momentum tensor ([2]) , respectively. In this approximation we have the following expressions for the total mass density 
and the pressure: 

Pt « fp, Pt « fp. 

Substituting these expressions into (|44[) and taking into account that p = p c 9 n and p = ap c 9 n+1 , we have from (|4"4"|) : 



dp„ 



nd9 



df 



(l + a6) f{i + aey 



This equation differs from the case without a scalar field by the presence of the second term on the right-hand side 
containing the coupling function /. In the absence of the non-minimal coupling, i.e. when / = 1, this term vanishes. 
Integrating, we find 



Pg = Pgc 



'{1 + <t)0~ 


n 


\ [ df 


exp 


1 + <j9 _ 


[J f(l+*0)_ 



(45) 



where p nr is the value of the gas density p g at the center of the configuration. This constant may be evaluated as 
follows [28|]: near the boundary of the configuration, where 9 -C 1, equation (|4"5|) becomes approximately 



Pg 



where £ a is the point starting from which the approximation 9 <C 1 becomes valid. This equation can be rewritten in 
differential form as 



d_ 



In 



Pg 



Pgc(l 



d 



(46) 



Since / = f(<fi(0) 1S a function of £, the term on the right-hand side can be evaluated as follows 



d_ 



*df_ 



d 



(df/dQdA 



= ^7 (In/). 



Taking this expression into account, equation f|46[) takes the form 

Pg tt Pgc {l + a) n f6 n , 

Near the boundary the internal energy density is small compared to the rest-mass energy density, so that pt 
Comparing the above equation with the expression p t ~ f p c 9 n , we obtain 



Pgc 



(l + d)"' 



and equation (|4"5"j) becomes 
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Using expression (|47|) . the corresponding proper energy Eq 9 of the gas may be defined as the integral of p g with 
respect to proper volume dV = knr 2 e x l 2 dr for the metric (|3]) as follows 

E 0g = M 0g = 4tt / Pg e x ' 2 r 2 dr. (48) 



The quantity Mq 9 is proportional to the total number of particles N in the configuration, Mq 9 = N/A, where A is 
Avogadro's number. Using the dimensionless variables (fTS]) and expressions (|30l) . (|47p . we obtain the proper energy 
of the gas l|48p in units of the total energy E = M in the following form 



'{it [/(i + ^r^/}e 2 ^ 

E ~ Jo [l + <j9} n [l-2<j{n+l)v/Z} 1/2 



E 0g _ 1 V ex P 



(49) 



This expression differs from the corresponding expression in [28j through the presence of an extra factor exp{...} under 
the integral which comes from the non-minimal coupling in the system. Stability of the configuration can tested via 
the sign of the expression [2^, l32j 



Binding Energy = E 0g - E = Eq^ 

E E ~ E 



1. (50) 



The necessary (but not a sufficient) condition for stability of the system that Eo g > E i.e. that the total system 
energy E is less that the energy of the non-gravitationally interacting gas, E$ g , making E the energetically preferred 
state. The condition in ((511)) amounts to requiring that the Binding Energy defined in this equation be positive. In 
the paper (34j |. the question of stability of adiabatic polytropic configurations was considered. It was shown that such 
configurations may have both positive and negative binding energies depending on the value of the polytropic index 
n, and the parameter a (configurations with 1 < n < 3 and 0.01 < a < 100 were considered). 

When a non-minimally coupled scalar field is included, we considered a very narrow range of these parameters, 
(7 = 0.2 and n = 1.0, 1.5, and also only one choice of the coupling function / given by the form (f2"2"| . In this case the 
numerical calculations indicate that Eo g /E <C 1, i.e. we have the negative binding energy, and correspondingly the 
configurations under consideration are unstable. The same values of a and n without a scalar field give configurations 



}■ 


\m)] 




[ /(0) J 



with the positive binding energy [34|. Obviously, the difference between the present results and those of [34| is 
connected with the presence of the extra factor exp{...} under the integral in the expression (|49[) whose value, as 
numerical calculations indicate, is much less than unity for the parameters a, n used in the paper, and for the 
coupling function / of the form (f2"2"|) . 

The behavior of exp{...} in (|49p can be approximately estimated as follows: Since we are looking only for solutions 
with < 6 < 1, the value of (1 + a6) may be approximated as giving a finite contribution to the value of integral 

J^ 1 [/(l + a 6)} 1 df in the form of some constant factor a > 0, i.e. J^ 1 [/(l + a 6)} 1 df ~ a J din/. Then one can 
sec that 

expj^ [/(l + aOT'df 

When the coupling function / is taken in the form (|2"2")l . we have the finite value of /(0), and since </>(£i) —> then 
— > as well, thereby suppressing the value of E 0g /E in (|4^|) . In the non-relativistic case, where a — > 0, we have 
a — y 1 and the situation with the instability is the same. 
Possible ways of resolving this problem are: 

(i) Consider configurations for which </>(£i) 7^ 0, and the function /(£) varies slowly in the range < £ < £1 to 
provide a large value of the above exp{...}. In this case, one can expect that the configuration will have a "tail" of 
a scalar field outside the fluid at £ > £1 (a similar thing happens for configurations considered in [f|). In this case, 
it will be necessary to perform a stability analysis of both internal solutions (where the fluid and the scalar field are 
involved) and external solutions (where only the scalar field is involved). 

(ii) Since the solutions and their stability are sensitive to the nature of the scalar field-fluid coupling, which in this 
paper was take to have the form given in (|22[) . one could look for stable solutions by examining different forms of this 
coupling; for example one could try a non-polynomial coupling of the form / = e^. These stability studies will be the 
focus of future work. 
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V. NON-RELATIVISTIC CASE: AN ANALYTICAL SOLUTION FOR A MASSLESS SCALAR FIELD 



The numerical results obtained in the previous sections demonstrate the possibility of obtaining the regular solutions 
both in relativistic and non-relativistic cases when the coupling function / is chosen in the form of (JUJ). This choice 
is not the only one possible. It will be shown latter in this section that, by choosing the special form of the function 
/, it is possible to find an analytical solution in a particular non-relativistic case when the polytropic index n = 0, 
and the scalar field is taken to be massless. This case corresponds to an incompressible fluid with a constant mass 
density, p = p c = const, and a spatial varying pressure, p ^ const. 

To begin we rewrite the system of equations (j2~4"]) -([2l) |) for the massless case and in the non-relativistic limit (a — > 0) 
in the form 



d? 



,d6_ 
di 

dv 

d£ 

2 



2 Vd£ 



2\d£ 



(51) 
(52) 
(53) 



We have kept the a-containing-tcrm on the right-hand side of equation (|53|) since it is does not necessarily small for 
arbitrary f. Differentiating the first equation of the system, and substituting the second and the third equations into 
this expression, gives 



(54) 



The last term of this equation contains the factor df/d£ = (fi'df /d<f>. For / = 1, equation (|54[) reduces to the Lane- 
Emdcn equation (|38[) which has an analytical solution for the case of an incompressible fluid with a constant mass 
density p = p c = const, p ^ const, n = - see (32j . equation (10.3.14). This solution is 



i--e 

6^ 



(55) 



Equation (|54[) can also be integrated analytically by choosing the function / to have the following power-law form 



/ = h + -a 



(56) 



where fo, /3 and m are arbitrary parameters. Using this function and taking the parameter m to be of order 1 so that 
the product a8 <?C |l/m|, we have from (fM]) : 



LA. 



,d6 



fo + -r 



When n = 0, this equation can be integrated giving the following regular solution 



fo 
6 



e 



m(m + 2)(m + 3) 



m+2 



(57) 



Here 9q is an integration constant. A second integration constant was set equal to zero to make the solution regular at 
£ = 0. We take 6q — 1 so that in the absence of the scalar field (i.e. when fo = l and j3 = 0) equation (|57|) will reduce 
to equation (|55p . This solution is regular everywhere, including the point £ = 0, when m > —2. The degenerate case 
m = corresponds to the classical solution (JSSJ) which can be obtained from (f5"T|) by setting fo = l and = 0. This 
choice that corresponds to turning equation (|54[) into the Lane-Emden equation (|38l) . 

We now analyze the behavior of the solution as £ — > 0. To begin we calculate the first and second derivatives of 8: 



dl 



ft 



m(m + 3) 



rn+l 



~d~e 



fo 
3 



g(m + l) 
m(m + 3) 



For a spherically symmetric solution with maximal mass density as £ — > 0, two conditions must be satisfied: (1) the 
first derivative must be equal to zero; (2) the second derivative must be negative. The former condition gives m > — 1. 
Given this the behavior of the second derivative is: 
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if (-1 < m < 0): 



-oo for (3 > 0; 
-oo for /3 < 0; 



if m > 0: 6»" -> -/ /3 for any /3. 

One can see that for (— 1 < m < 0) the second condition is satisfied only if j3 < 0. 

On the other hand, the regularity of the solution assumes the presence of a point on the axis, £, where 8 = 0. For 
the classical solution (|55|). this point is £i|e=o = \/6. In the case of configurations with a scalar field, the location of 
the point £i|e=o w ^ depend on the factors in front of the £ m + 2 term in ([57)) . i.e. it will be determined by the values 
of P and the denominator m(m + 2)(m + 3). In the range ( — 1 < m < 0) only negative /3 are allowed if one wants 
9" < as £ — > 0, which is one of conditions that the solution ([57)1 be regular. In the range m > the regularity of 
the solution is guaranteed if /? > 0, since otherwise the solution grow as a power law of £. 

The above analysis shows that the inclusion of a scalar field decreases the size of configurations, as compared to 
the classical case ([55)). for all acceptable values of the parameter m. In this sense the behavior of the solutions with a 
masslcss scalar field and the coupling function in the form of (|56[) differs from the behavior of the solutions obtained 
in section HV B I when the function / was chosen to have the form ([2"2")) . In the latter case, the size of the configurations 
for given values of n and a depends considerably on the parameters A, /3, and <fio, and it can be larger or smaller than 
the size of a configuration without a scalar field. 

Substituting the expression for 9 from (|57[) into the scalar field equation ([53")) with n = 0, we find the following 
analytical solution for the scalar field <f> in the interior of the configuration 



_ 

6 



m 2+m m 2 



m(m + 2) 2 (m + 3) (2m + 3) (m + 3)(m + 4) (m + l)(m + 2) 



C n+ • (58) 



One can see from this expression that the behavior of the scalar field is defined in large measure by the factor (3a in 
front of the square brackets. Since in the non-relativistic limit a — >• 0, this factor will differ appreciably from zero 
only for large j3. For (3 not large one can expect that the field will be practically constant and equal to the central 
value cf>o up to the boundary of the fluid at £ = £i where 9 = 0. The general solution for <j> for all £ is obtained by 
matching of the internal solution (|58j) with an external solution of equation (|53j) where the right-hand side is equal 
to zero. The exterior equation has a solution of the form 

<t>ext = Ci + C 2 /£, 

where C± , Ci are integration constants which are determined by matching of the external solution <f> ex t with the interior 
solution 4> in solution at the boundary of the fluid at £ = £i. One can see from this expression that the corresponding 
mass density, which is proportional to 4>'^ xt , tends asymptotically to zero. 



VI. CONCLUSION 



In this article wc studied gravitating, spherically symmetric, star-like, configurations with a matter source consisting 
of a normal (i.e. non-ghost, non-phantom) scalar field plus a perfect isotropic fluid. The motivation for studying such 
a model is that scalar fields have found broad use in various cosmological models as well as astrophysical models. It is 
natural to postulate that, if scalar fields do really play a role in the Universe that these scalar fields might play a role 
in the structure of compact objects such as Main Sequence stars or neutron stars. Proceeding from this assumption, 
we studied a model of a star-like configuration supported by a scalar field non-minimally coupled to ordinary matter 
in the form of a perfect fluid. As an example, we considered the case when the scalar field had a quadratic mass 
term and a quartic self-interaction giving a scalar field potential of the form (|23[) . The coupling between the scalar 
field and the perfect fluid, /, was taken to have the form ([2"2")) . For this model we studied solutions both in the 
rclativistic and non-relativistic limit. Our investigation showed that the existence of regular solutions is possible 
exactly because of the presence of the non-minimal coupling between the scalar field and fluid. In the absence of such 
coupling, the potential ([2"3")) gives only singular solutions From the mathematical point of view, the existence of 
regular solutions in our model is possible because of the appearance in scalar field equation (|26p of an effective mass 
term m e // = (fi 2 — (3a9 n+1 ) whose sign depends both on the behavior of the fluid density 9 and the values of the 
parameters /3 and a. It was shown that when m e ff < 0, there were regular solutions with finite masses and sizes. 
In some sense the solutions presented above are a cross between the interacting scalar field solutions of [33| and the 
interacting real scalar field solutions of [3]. In the interior region, < £ < £i, we have a real scalar field and a fluid 
which interacts with the scalar field. Due to the chameleon-like behavior of the scalar field with respect to the fluid 
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(from figure [3] one can see that in the interior region the field mimics the fluid) one might think of the two degrees of 
freedom associated with the scalar field and the fluid as being equivalent to the two degrees of freedom of a complex 
scalar field. While there is certainly some validity in this analogy (both the present solutions and those in [33[ are 
non-singular) one must avoid pushing the analogy too far since the fluid vanishes exactly at some point (i.e. 9 = at 
£ = £i) while the complex scalar field of [33[ and the real scalar of the present solution only go to zero asymptotically. 
Thus after reaching the point £ = £i our solutions become those of [lj] . 

Our results are interesting since they show tha t by adding such a fluid one can find regular solutions, thus evading 
some gravitational version of Derrick's theorem [23l . |24| which prohibits regular solutions for the system of gravity 
plus a normal scalar field if the potential V(ip) > 0. Thus the addition of the fluid was crucial to the existence of these 
solutions. The original, non-gravitational version of Derrick's theorem [U [22j assumes not merely the regularity of 
solutions, but also their stability. In this paper we have performed a preliminary stability analysis based on energy 
considerations. In this case we compared the total energy of the system (including the internal and gravitational 
energies) and the rest energy of the gas particles. For the values of the polytropic index n and the parameter a, 
and the specific coupling function / in the form of (|22|). this energy approach showed that the solutions studied in 
this paper are unstable. This was because the binding energy of the system, which was equal to the difference of 
the rest energy and the total energy, was negative. But, as it was shown in [34| . even a positive binding energy does 
not guarantee stability of a system allowing its transition to an energetically more advantageous state with the same 
polytropic index n, but having another parameter a. Obviously, such process is accompanied by an ejection of the 
excess of energy. In our case, the possibility is not excluded that there exist regular solutions with positive binding 
energy, for other values of n and a and other parameters of the scalar field /3, A, <pg. There is an additional possibility 
for finding stable solutions which was described at the end of section HV CI - since the solutions and their stability are 
sensitive to the form of the scalar field-fluid coupling (which in the present work had the form (|22[) ) one could look 
for stable solutions by changing the form of the scalar field-fluid coupling. For example, one could try non-polynomial 
couplings of the form /(</>) = e*". Finally, the energy approach to stability which was used in this paper, can and 
should be supplemented by a study of dynamical stability along the lines of the linear stability analysis performed in 
[3 EM or by using the catastrophe theory method suggested in [35| . In future work we plan to perform a stability 
analysis of the chameleon star solutions investigated here using both these approaches. 

The spherically symmetric solution that we found for the system considered in this paper were called chameleon 
stars in analogy with chameleon cosmological models [l(| because the characteristics of the scalar field (e.g. its mass) 
strongly depend the other fields and fluids in its environment. We characterized the behavior of these chameleon 
star configurations for both rclativistic and non-relativistic cases. For the relativistic chameleon stars we found that 
in general the mass-energy density tended to be more concentrated toward the center of the star at £ = 0. Despite 
this greater concentration of mass near the center of the chameleon star the total mass of the chameleon star was 
lower than the corresponding relativistic stars without a scalar field. This was the result of the outer regions of the 
chameleon star being much less dense than a corresponding non-chameleon, relativistic star. Similar comments apply 
to the non-relativistic chameleon stars. For the non-relativistic case we were able to find an analytical solution for 
the case when the scalar field was massless. This analytical solution had similar behavior to the numerically obtained 
solutions from section HVBI 

We briefly give some speculations about possible physical applications of the present solutions. Scalar fields are 
thought to play a role in cosmological dynamics (dark energy) and in dynamics at the galactic scale (dark matter). 
The present proposal is that a scalar field could play some role in astrophysics at the stellar scale in the formation of 
the chameleon star configurations discussed here. The fluid in our model would be provided by the star. For example 
one could consider neutron stars which contain a significant amount of the scalar field. From the discussion above 
our chameleon stars - both relativistic and non-relativistic - would be less massive for a given radius despite having 
a higher concentration of mass near the center. Thus one might look for neutron stars which have a larger radius 
for a given mass than would be expected for a normal neutron star. In this sense chameleon stars would have the 
opposite behavior from hypothetical quark stars which have a larger mass for a given radius. Another possibility is 
that some living stars (i.e. stars which are still fusing elements and are on the Main Sequence) might have trapped 
some significant amount of scalar field in their interior. Such living, chameleon stars would tend to have a larger 
density near their center thus increasing the rate at which they fused elements. The consequence of this is that such 
stars would be hotter and live a shorter than expected life span for their mass. Again a given chameleon star would 
have less overall mass than a non-chameleon star but would be nevertheless have a higher, interior temperature. The 
physical applications of the chameleon star model suggested above would of course require the existence of stable 
solutions. A search for such stable solutions is a goal of future studies. 
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Appendix A: Derivation of the energy-momentum tensor 

In the case of the isentropic quasipotential flow, the Lagrangian of the continuous medium has the form [29l |3C 



L m = p = W - e, 

where W = (e + p) is the heat function, and e is the energy density. Introducing the quasipotential 

dS 

cSi = c—— = Wui, 
ox 1 

where c is the velocity of light and S is the action for the matter, we have 

i w 2 i W 2 
S[S — —iiiU = — , 

c c z 

from which it follows 



(A.l) 



W =cy/g ik SiS k . 



The Lagrangian (|A. 1|) now takes the form 



L m = cy / g ik S l S k - e. 

The total action (including gravity plus all matter sources) in curvilinear coordinates is [If 



S = - [ L^f—gd^x, 
c J 



(A.2) 



(A.3) 



where the Lagrange density, L, above is that given in equation (JTJ. By varying the action (|A.3j) with respect to 
the metric, gik, one can obtain the Einstein equations and the energy-momentum tensor. In this way the energy- 
momentum tensor is given by the expression 



2 V 9 lk ~ dg lk 8x l d 9£l ' 



where the Lagrangian, £, contains only the matter components of the Lagrangian ([T]) namely 

£ = \g lk d lV d W ~ V(cp) + f(<p)L m . 



(A.4) 



Since v / ~ff£ does not depend on j?-r, the last term in (|A.4[) equals zero. Then we have 



1 /— r . = d^—gZ ^d^g 
2 V 9 lk i)g- dg* k 



-d£ 
dg 



(A.5) 



Taking into account that dg = gg tk dgik = —ggikdg lk and also using the expressions (| A. 1 1) and (|A.2|) . we find 

8^—g 1 



and 



- d£ 



1 cSiS k , , 1„ „ 



-g (fWu.Uk + diipd k ip) . 
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Substituting these expressions into (|A.5|) . we finally have: 

Tik =/[(£ + p)iHUk - g ik p] + diipd k (fi - 9ik 



-d^if - V(cp) 



(A.6) 
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